Two-photon decays of tt^, rj and 77' 



B. Borasoy\ R. Nifiler^ 

Physik Department 
Technische Universitat Miinchen 
D-85747 Garching, Germany 



Abstract 

We investigate the decays of vr'', 77 and r]' into two photons in an effective U{3) chiral 
Lagrangian approach without employing large Nc arguments. Tree level and one-loop 
contributions from the anomalous Wess-Zumino-Witten Lagrangian are calculated and 
the importance of ij-rj' mixing is discussed. Unitarity corrections beyond one-loop play an 
important role for the decays with off-shell photons and are included by employing a cou- 
pled channel Bethe-Salpeter equation which satisfies unitarity constraints and generates 
vector-mesons from composed states of two pseudoscalar mesons. 
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1 Introduction 



The two-photon decays 7r°, ry, r^' — > 77 are phenomenological manifestations of the chiral anomaly 
of QCD and can provide important information on the chiral symmetry of the strong interac- 
tions. In order to study the phenomenological implications of the chiral anomaly, it is convenient 
to employ a chiral effective Lagrangian, since at low energies the effective degrees of freedom 
are colorless hadrons rather than quarks and gluons. The effective Lagrangian contains a piece 
which reproduces the anomalous behavior of the effective action under chiral transformations. 
Such a Lagrangian was systematically constructed by Wess and Zumino [1] by directly integrat- 
ing the anomalous Ward identities and later Witten provided a representation of the anomaly 
which illustrated the topological content of the theory [2]. 

There are many anomalous processes which can be calculated from the Wess-Zumino- Witten 
(WZW) Lagrangian, such as 7r° — > 77, r) — > 77, 7 — > Stt, r) — > tt'^tt'^ etc. Originally, the 
WZW Lagrangian was formulated for the eight Goldstone bosons (vr, K, rj) which form an octet 
under chiral S'f/(3) transformations, but it can be easily extended to include the 77', the singlet 
counterpart of the Goldstone boson octet [3, 4, 5, 6]. Although the rj' is not a Goldstone boson 
due to the axial U{1) anomaly of the strong interactions, it combines with the Goldstone bosons 
to a nonet at the level of the effective theory. This allows the phenomenological investigation 
of r]' decays, such as r]' —>■ 77 and rj' — > 7r"'"7r~7. 

The present work deals with the two-photon decays of tt", t] and rj' with one, both or none of 
the photons being off-shell, 7i'^,ri,ri' — >• 7(*)7(*). The decays P ^ 'yy (P = 71^,7], rj') with both 
photons on-shell are dominated by the WZW Lagrangian. They have been calculated up to one- 
loop order within chiral perturbation theory (ChPT) and were shown not to receive non-analytic 
corrections from the one-loop diagrams which were compensated by the chiral corrections of 
the pseudoscalar decay constants [3, 4]. This is also in agreement with the complete one-loop 
renormalization of the anomalous Lagrangian [7, 8, 9], where it was shown by using heat kernel 
techniques that in dimensional regularization one-loop diagrams for P — > 77 do not lead to 
divergences that must be renormalized by appropriate counterterms of higher chiral order. It 
was furthermore argued [3] that a consistent picture of 77-77' mixing emerged from the two- photon 
decays with one mixing angle of approximately —20°. 

More recently, a two-mixing angle scheme has been proposed by Kaiser and Leutwyler 
[5, 10, 11] for the calculation of the pseudoscalar decay constants in large N^. chiral perturbation 
theory. The two angle scenario has been adopted in a phenomenological analysis of the two- 
photon decay widths of the r) and rj', the 777 and 7/'7 transition form factors, radiative J/^ 
decays, as well as of the decay constants of the pseudoscalar mesons [12, 13]. The authors 
observe that within their phenomenological approach the assumption of one mixing angle is 
not in agreement with experiment whereas the two-mixing angle scheme leads to a very good 
description of the data. As pointed out in these investigations the analysis with two different 
mixing angles leads to a more coherent picture than the canonical treatment with a single angle. 
In particular, the calculation of the pseudoscalar decay constants within the framework of large 
Nc chiral perturbation theory requires two different mixing angles [10] . (A similar investigation 
was performed in [14] but with a different parameterization.) 

In addition, it was shown in [15] that even at leading order 77-77' mixing does not obey 
the usually assumed one- mixing- angle scheme, if large Nc counting rules are not employed. 
One purpose of this work is to critically re-examine the contributions from 77-77' mixing to the 
two-photon decays. 
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For the decays P — > 77* contributions from higher chiral orders will become increasingly 
more important for larger photon virtualities. In this case one-loop contributions turn out to be 
divergent and must be renormalized by appropriate counterterms [4, 7, 8]. The remaining finite 
parts of the relevant counterterms can be estimated by vector-meson exchange contributions 
[8]. In [8] the vector- mesons were included explicitly and integrating them out from the effective 
theory produced contact interactions. An underlying assumption of this approach is that the 
masses of the vector-mesons are much larger than the involved momenta. 

Alternatively, the chiral Lagrangian for the Goldstone bosons is able to reproduce a number 
of meson resonances such as the /o(980) and the p(770), when combined with non-perturbative 
Bethe-Salpeter approaches which are employed in such a way so that they ensure unitarity 
[16, 17]. Within these approaches effective coupled channel potentials are derived from the 
chiral meson Lagrangian and iterated in Bethe-Salpeter equations (BSEs). The BSE generates 
dynamically quasi-bound states of the mesons and accounts for the exchange of resonances 
without including them explicitly. The usefulness of this approach lies in the fact that from a 
small set of parameters a large variety of data can be explained and no additional assumptions 
need to be made on the couplings of the resonances. In [18] the approach has been extended 
to include the rj' and successfully applied to the hadronic decays of rj and rj' in [19]. The 
main purpose of the present work is to embed the coupled channel approach in the two-photon 
decays P — > 7(*)7(*) and to investigate the importance of vector-mesons for these decays from 
a different perspective. This will also allow us to give predictions for the double Dalitz decays 
with two off-shell photons which have not yet been studied experimentally. It may further 
clarify the question whether double vector-meson dominance holds, which is also an important 
issue for the anomalous magnetic moment of the muon and kaon decays [20]. 

This work is organized as follows. Section 2 introduces the notation and the tree level 
contributions to the decays which arise from the WZW Lagrangian. One-loop contributions 
and mixing effects are described and discussed in Sec. 3. The application of the coupled channel 
analysis to the decays is outlined in Sec. 4 and the numerical results are presented in Sec. 5. 
We conclude with a summary of our findings. 



2 Decays at tree level 

The effective Lagrangian for the pseudoscalar meson nonet (tt, rjs, rjo) reads up to second 
order in the derivative expansion [5, 11, 21, 22] 

£(0+2) ^ ^ ^^^^^^t£,Mc/) + v^iU^'x + X'U) + iVsiU^'x - X^U) 

+ V^{U^D''U){U^D^U) (1) 

where wc have presented only the terms relevant for the present work. The unitary matrix U is 
a 3 X 3 matrix containing the Goldstone boson octet (tt^, tt^, K, rjs) and the rjo. Its dependence 
on 7r°, 778 and rjo is given by 

U = exp(diag(l, -1, 0) • zttV/ + diag(l, 1, -2) • irjs/Vsf + iV2rio/V3f + ...). (2) 

The expression (. . .) denotes the trace in flavor space, / is the pion decay constant in the chiral 
limit and the quark mass matrix M — diag(mu, m^, m^) enters in the combination x — 2SA1 
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with B = — (0|gg|0)//^ being the order parameter of the spontaneous symmetry violation. The 
covariant derivative is defined by 



D^U^d^U -i{v^ + d^)U + iU{v^-d^). (3) 

The dependence of the effective Lagrangian on the running scale of QCD due to the anomalous 
dimension of the singlet axial current A° = |?7/i75? is absorbed into the factor y/X in the 
axial- vector connection 

at, = a^ + — — — (a^). (4) 

which is the scale independent combination of the octet and singlet parts of the external axial- 
vector field ttfj,, cf. [22] for details. Due to its scale dependence, \/X cannot be determined from 
experiment, and all quantities involving it are unphysical. 

The coefficients Vi are functions of rjo, 1^(770//), and can be expanded in terms of this 
variable. At a given order of derivatives of the meson fields U and insertions of the quark mass 
matrix one obtains an infinite string of increasing powers of rjQ with couplings which are 
not fixed by chiral symmetry. Parity conservation implies that the Vi are all even functions 
of r]o except V3, which is odd, and Vi{0) — ¥2(0) = Vi{0) — 31^4(0) = ^/^ gives the correct 
normalization for the quadratic terms of the mesons. The potentials Vi are expanded in the 
singlet field 770 



Vi 
V3 



i2 ^4 



f]=«r+4^f + «f^ + ... for i = 0,l,2.4 



/ 



,a)^ + 43)^ + ... (5) 



with expansion coefficients fp"* to be determined phenomenologically. 

The Lagrangian contains only terms of natural parity. The photonic decays tt — > 

77,77 — > 77,77' — > 77, on the other hand, are not covered by since they arise from the 

unnatural parity part of the effective Lagrangian which collects the terms that are proportional 
to the tensor e^^ap- Within the effective theory the chiral anomalies of the underlying QCD 
Lagrangian are accounted for by the WZW term [1, 2, 5] ^ 

where 

W{U, v) = {EUUvUv + Evdv + Edvv - iY?v) - {U ^ C/^) (7) 

with E = U^dU and we set Nc — S for the number of colors. We have displayed only the pieces 
of the Lagrangian relevant for the present work and adopted the differential form notation of 
[5] 

V = dx'^v^, d = dx^d^ (8) 

^Note that for our purposes we can safely set the singlet axial vector field (a^) and the derivative of the 
QCD vacuum angle, 8^9, to zero in Swzw which enables us to work with the renormalization group invariant 
form of the anomaly. 
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fc, e 




fc', e' 

Figure 1: Tree diagram of the decay P — > ^(*)^(*), where k^'^ and e^'^ denote the momenta and 
polarizations of the photons, respectively. 

with the Grassmann variables dx^ which yield the volume element dx^dx'^dx"dx^ = e^^°^^d'^x. 
The operation [U <-> U^) indicates the interchange of U and . The first integral on the 
right-hand side in Eq. (6) is taken over the five-dimensional manifold M5 which is the direct 
product of the Minkowski space M4 and the finite interval < a;^ < 1, while the Grassmann 
algebra is supplemented by the fifth element dx^ . The matrix- valued fields U{x^x^) in the first 
integral are functions on this five-dimensional manifold and interpolate smoothly between the 
identity matrix 0) = 1 and the matrix U{x) = U{x, 1). The value of the first integral 
is independent on the particular choice of the interpolating functions U{x,x^), whereas the 
integration for the second term in Eq. (6) extends only over the Minkowski space M/i. The 
external vector field v = —eQA describes the coupling of the photon field A = dx^A^ to the 
mesons with Q = idiag(2, —1, —1) being the charge matrix of the light quarks. 

At tree level only the terms quadratic in the vector fields v contribute and one arrives at 
the following pieces of the WZW Lagrangian 

d^xCy,^^ = ^—{UUUUUvUv + UUUvdv + UUUdvv 

-UdU^'UdvU^'v - UdU^'vdv - UdUUvv) + ... . (9) 

However, this is not the whole story, since there exist further terms at fourth chiral order which 
are gauge invariant and do not contribute to the chiral anomalies. A complete set of these 
terms has been given in [5], out of which the following three contact terms contribute 

d^xC^^ = Wi{UdvUUv) + W2{dvdv) + iW^idUdUUv + dUUUdv), (10) 

where the potentials VFi,2,3 are odd functions of the singlet field rjo. The first two terms con- 
tribute already at tree level, whereas the last term represents a vertex for a one-loop diagram 
and will be discussed in the next section. Following the steps of [22] it is straightforward to see 
that these terms are needed to absorb the QCD renormalization scale dependence via a set of 
redefinitions of the potentials Wj into the factor y/X of d^ in Eq. (4) . 

The decay at tree level is depicted in Figure 1 and the pertinent amplitudes are given by 

^(tree)^p ^ ^(*)^(*)) ^ k^e^k'^e^^^ j^a^'''^ (11) 



with 



47r2/ 
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k, e 
k', e' 

(b) 

Figure 2: One-loop diagrams contributing to P — > 'y^*^'y^*'>. In (b) the crossed diagram is not 
shown. 

where \ w^^ are the leading expansion coefficients of the potentials Wi and W2, respectively. 
It is important to note that in our scheme r)-r)' mixing is of second chiral order [15] and will 
modify the sub-leading order, i.e. the one-loop order of the decay amplitudes. This is in 
contradistinction to large ChPT where rj-r}' mixing contributes at leading order. 




3 One-loop contributions 

The one-loop diagrams for P ^(*)^(*) of order 0{p^) have anomalous vertices from the 
following pieces of the WZW Lagrangian 

d^xCy,^^ = — ^{UUUUUvUv + UUUvdv + UUUdvv - iUUUUUUUUUv) 

-{U ^U^) + ... . (13) 

After expanding U = exp(^0) in the meson fields cj) we arrive at the expression 

d^xU^^ = -:^^[{d(p [0, [(P,vdv]])- {d(P [0 , dv] [0, v]) - 2i{d(P d(P d(P v)') + . . . , (14) 

where we have only presented the relevant terms for the one-loop diagrams. The one-loop 
diagrams at order 0{p^) contributing to the decays P ^(*)^(*) are shown in Fig. 2. 
The tadpole contributions in Fig. 2a read 

A(^<^<^)^P ^ ^(*)^(*)) ^ e%e.k'^e'f,e^^'^^^^ (/?r^'^A. + /jf -^^'^Ax) (15) 
with being the finite parts of the tadpoles 

J i27rY P-ml + ie IGtt^ //^ ' y^-^) 

^ / finite 

and fi is the scale introduced in dimensional regularization. In the present work we are only 
concerned with the finite pieces of the diagrams and neglect the divergent portions throughout. 
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The coefficients p^^'"''-^''^ are given by 



n(tad),iT 
n(tad),K 



3' 
5 

3' 



(tad),- 



Pi,- 



(tad),K 



2 

1 

7!' 



(tad),iT 



^(tad),K ^ 



.2.„(1) 



(17) 



We have furthermore replaced the pion decay constant in the chiral hmit, /, by the expressions 
Fp which include the one-loop corrections. The corrections have been calculated in f/(3) ChPT 
without imposing large Nc counting rules and in infrared regularization [15]. We can employ 
the same formulae also for the present work by noting that only tadpoles contribute at one-loop 
order to the decay constants. In infrared regularization the rj' tadpole vanishes, whereas the 
tadpoles for the Goldstone boson octet remain unaltered. This implies that the rj' tadpole docs 
not contain any infrared physics and can be absorbed completely into the low-energy constants 
(LECs) of the effective Lagrangian. It is neither a function of the Goldstone boson masses 
nor of the external momenta, it is just a constant. We will therefore assume that rj' tadpole 
contributions have been compensated by renormalizing the LECs appropriately and will work 
with the renormalized values without indicating it explicitly. 

The expansions of the decay constants in terms of the Goldstone boson masses up to one- 
loop order are given here for completeness [15] 



F^ = f 1 + 4A 



3(0)- 



2mj^ + m: 



- + 



(0) 



P 



f 



l + 4/5f^^^^+4/3f^ 



P 



The LECs jS'f'^ originate from the natural parity part of the effective Lagrangian at fourth chiral 
order 

>C^'^ = Efc/3feOfc, (19) 
where the contributing fourth order operators are 



/ 



F-n'Q — / 



1 + 4 



2m]^ + ml 



^(3/3f -9/3i",^+3/3: 



3/2 



(0) 



?(0) 
18 



+3/?£^ + 9/?l?-3/?(?-|V6/?g)) 



(18) 



4 = -(C'^C^)(M), O5 
Oi7 = (C^)(C^)(M), 

046 = 2i(5^)(C^M), O47 

052 = -2(M)a'^(a^), O53 

and we made use of the following abbreviations 



-{C^){C,M), 

2i(a^)(C^)(M), 

2i{N)d^^{a,), 



(20) 



(21) 



The Pk are functions of the singlet field, Pkivo)^ ^.nd can be expanded in the same manner as 
the Vi{rio). Note that the also have divergent pieces, in order to compensate the divergences 
from the loops, but in the present work we are only concerned with the finite parts and use the 
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same notation for simplicity. The decay constant of the r( related to the singlet axial-vector 
current, F^/q, is defined by the matrix element 

(0|^g7M75g|?7')=iPM^r,'o. (22) 

Due to the anomalous dimension of the singlet axial current, F^'o depends on the running scale 
of QCD and its value cannot be determined phenomenologically. In order to account for the 
chiral corrections for /, we have therefore employed the scale invariant ratio F^i — -j^F^iQ in 
the expression for the rj' decay amplitude in Eq. (15). 

We replace the pion decay constant in the chiral limit, /, by the chirally corrected decay 
constants Fp also in the tree level result, Eq. (11). This replacement yields corrections at 
0{p^) which must be compensated. If the decay constants Fp in Eq. (18) are written as 
Fp — f{l -\- 5Fp/ p), it induces at sixth chiral order the corrections 



^(/)(P ^ y*)^W) = eXe.k'J,e^^-^^yr^5Fp. (23) 
The contributions from diagram 2b yield 

+(5(r^^''[I{ml;e) + I{ml-X)]). (24) 

The integral / is given by 

/(m^;/) = ^(^A + (m^ - ^)G^^(/) + ^(/ - 6m^)), (25) 
where G is the finite part of the scalar one-loop integral 



which reads 



167r2 



, mm m^ — m^ , m 
1 + In ^ + In — 

jjL'^ m 



2\/A~(F) ^^^^^j^ \/Amm(p2) 



p2 _|_ — 

^mrniv") = {{m - Tuf - p") {{m + TTlf - p") . (27) 

For the particular case with on-shell photons the integral / reduces to the simple expression 
/(m^; 0) = A. The coefficients are given by 

i(uni),TV n(uni),TT 1 r)(uni),'K /2 ic 2 (1) 



V V 3 (2g) 

(uni),K _ -. a(um),K _ a(uni),K _ ic^2^„(l) 
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3.1 Wave-function renormalization, 77-77' mixing, and counterterms 

There are further chiral corrections at 0{p^) which arise from wave- function renormahzation 
and r]-ri' mixing. The relation between the original SU (3) fields (vf^, rjs, i]o) and the renormalized 
states (tt",?],?]') without employing large Nc counting rules has already been derived in [15]. It 
is given in terms of the 3x3 matrix 1 + i?^^), (7f°, rjs, 770)^ = (1 + i?(^))(7r°, 77, f]')^, with 

< = (-4[2m^ + mllpP - 4mJ/3f + iA;,)/F,^ 

< = ^« - rn^) (2^^/^5,18 - 4'')/Fiml 

- + mlmT + - 9/3i? + ^P[^)/Fi, (29) 

and all remaining entries vanish ^. In Eq. (29) we have used the abbreviations 



° 2 


2.f 






P ' 














3g{0) 



(30) 

Within this scenario, rj-r]' mixing contributes at next-to- leading order which is in contradistinc- 
tion to large ChPT [23], where it is a leading order effect. Inserting these relations into the 
tree level result Eq. (11) one obtains the following corrections at sixth chiral order 



^(^)(P 7(*)7(*)) = e%e.k'J,e^'^'^^^^/3^'^ (31) 

with the coefficients 



(2) 
■Ori 



a{Z) _ (tree)T^i2) 

/5f =«r^<UaJ-X'- (32) 

Summing all the contributions at order 0{p^) we verify that for the vr"- and r^-decays into two 
on-shell photons the dependence on the regularization scale fi cancels out, i.e. the amplitudes 
are finite and do not need to be renormalized by counterterms of the unnatural parity part of 
the Lagrangian [3, 4, 7, 8]. For the process rj' — > 77, to the contrary, the w^^ term of the 
unnatural parity Lagrangian at fourth chiral order induces a non-vanishing fj, dependence via a 
loop. Also, for the decays with one or both photons off-shell the divergent parts do not cancel 
for any of the decaying mesons and again a residual scale dependence remains which can be 
compensated by introducing counterterms at sixth order. The explicit renormalization of the 

^Note that we work in the isospin Umit so that the 7r° field does not undergo mixing with the rjs-rio system. 
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amplitudes which are induced by the WZW term has been accomphshed in [7, 8] and is beyond 
the scope of the present investigation, where we restrict ourselves to the finite pieces of the 
integrals and counterterms. 

In the SU{3) framework without an explicit rj', the Lagrangian of unnatural parity at 0{p^) 
has been constructed in [24]. The relevant terms contributing via tree diagrams to the decays 
are presented in App. A, where we also show an additional set of counterterms which arises 
due to the extension to the U{3) framework. These terms can be reduced to the following 
structures (restricting ourselves to the finite pieces of the coupling constants) 



,(6) 



dv □ dv) 



w. 



(1) 16 



rjQ (dv □ dv) 



+ ^^^^dvdv) +wf^ ^{cl)x){dvdv) 
+ w^^^ jTjoixdvdv) +w^e^ ^ rjo{x) {dv dv) , 



(33) 



where we employed the LECs of the original counterterms in App. A. Note that the first two 
terms only contribute to the decays with off-shell photons. The contribution of £^^^ to the 
decays reads 



with the coefficients 



(34) 



Pr,' — 



64 



(0) 



3 
64 



3v/3 
12 



TT 



-TT 



+4(2m^ + m^)(;/^«;f 



4, 2 



32 (m^ 



2ml) { 



ml)w\ 



w 



(0) 



w 



w, 



(35) 



3.2 Numerical results at one- loop order 

Before implementing the unitarity corrections beyond one-loop within the Bethe Salpeter ap- 
proach, we would like to extract numerical results from the one-loop expressions for the decays 
into two physical photons, P 77. Combining all the contributions of the preceding sections 
we arrive at 

^(l-loop) _ J^(tree) _^ J^(tad) _^ J^(uni) _^ j^{f) _^ j^{Z) _^ ^(ct)_ ^3g^ 

The decay width F is given by 
with a — e^/47r and 

= a^*'--^ ( 1 + + '^^.^ + 2/?^^"^^'")^ + (/?£'^'^'^ + 2/?^^"*^'^)^ + I3f^ + 13^^'\ 

\ Fp J Fp Fp 

(38) 
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For comparison with previous work [3, 4] we also present the exphcit form of the (5p~^""^^ 
coefficients 



q(I-Ioop) 
p(l-loop) 



1 - 


2567r^ 
3 


2 -(0) 


1 




1 32 


7! 





1 - 



(m^ — m 



'ttJ o 



2 



+ [(4m^ - 7m^)^f + 24(mi - m^)^f ]} , 



(l-loop) 



F2 



+ 



512 



167r^(^?) - 2^.«)(A. + A^) + W| ^(m^ - m^)(2/35,i8 - 4^ 



-n 



rriQ 

« + 2m^)(J|^f + 4^)) + 3(2m^ + m^)(w|«;f + 4^^; 



(39) 



For the pion decay constant we take the value F^o = = 92.4 MeV, while for we employ 
1.3F^ which follows from an analysis within the framework of conventional SU{Z) ChPT 
that is similar to our expression at the order we arc working [21]. Since both the values of 
the contact terms w^l\ w^2^ and wl^\ which contribute to the rj' decay as well as to the rj 
decay due to rj-r]' mixing, and the values of the 0{p^) couplings wf^\...,wf\id''^\w^Q^ are 
unknown, we first set them by hand to zero in order to see, if a reasonable fit to all three decay 
widths is possible without them. The difference between both mixing angles for 77-77' mixing is 



proportional to the parameter combination (3^^ + [15] which is dominated by f3^' , since 
/3i8^ represents an OZI violating correction. We will thus neglect (3^^ and borrow the value 
for (3^^ from conventional SU{3) ChPT, (3^^ — 1.4 • 10~^ [25]. Furthermore, the parameter 

combination + SP^^ — — ^J^P^^ represents a, 1/Nc suppressed contribution to F^' in 

Eq. (18) with respect to P^^ + ^P^^^ — Sp'ff + p'^ . The value of the latter combination has been 

estimated in [15] to be roughly 0.15 x 10~^ and therefore pf^ +SpP —P^^s — \J\p^2 expected 

to yield a tiny correction to P^!~^°°^\ so that for our purposes it can be set to zero. The only 

remaining parameters are then F^/ and v^^ which we fit to the decay widths T^o — 7.74 eV, 
= 0.46 keV, and F^/ = 4.28 keV. These are the central experimental numbers quoted by the 
Particle Data Group [26]. The fit yields F^/ = 1.29F^ and v^^^ = 1.15F2/4. The value for v^^^ is 
in agreement with the one derived from an analysis of the pseudoscalar meson masses and decay 
constants [15]. The extracted value for F^/, on the other hand, is slightly larger than the values 
deduced within the one-mixing angle scheme, F^i pa I.IF^ [3, 4]. This would indicate that in 
the present approach some of the OZI violating corrections for F^/ in Eq. (18) are comparable 
in size with the leading contribution P^\ If, on the other hand, one fixes F^/ ^ I.IF^, a fit with 
only v^^ as a free parameter is no longer possible even within the experimental error ranges. 
We also compare our results with the one-mixing angle scheme by putting P^^ to zero. By 



?(o) 



fitting F^/ and v^^ to the decay widths we obtain F^/ 



1.15F^ and 



1.15F2/4 which 
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would correspond to a mixing angle of about —20° in nice agreement with [3, 4]. If we include 
the contact terms w^^ and wl^\ we have the freedom to choose F^' = 1.1 while still 
being able to match the experimental data. Changing the LECs w^]^\ w^^ and w^^ within a 
reasonable range of —3.0 x 10~^ ... 3.0 x 10~^ that is motivated by large Nc considerations and 
comparison with the coefficients of the WZW term we find a variation of the mixing parameter 
of V2 = (1.25 ± 0.18)F^/4 which is commensurable with a mixing angle of —21.3° ± 3.1° in 
the one-mixing angle scheme. The lack of knowledge of the exact values of w^^ and w^^ 
induces a 15 % uncertainty in the determination of the ?7-?7'-mixing angle. 

Since with the chosen parameters we are able to accommodate the experimental decay 
widths, while being in agreement with results from the conventional SU (3) framework, there is 
no indication that the unknown parameters of sixth chiral order wf'\ . . . , wf\w^^\wl^^ should 
have values significantly different from zero which would lead to sizeable contributions for the 
decay widths. Within the one-loop calculation presented here, we can thus neglect them. 

Note also that we work in the isospin limit of equal up- and down-quark masses. The effect 
of isospin breaking is only important for the decay of the 7r°. It has been discussed in [23, 27, 28] 
and yields a correction of about 5 % to the r^^o decay width. 



4 Unitarity corrections beyond one-loop 

In the preceding sections, we have identified the different types of contributions which arise 
from both the WZW and the unnatural parity p*^ Lagrangian at tree level and next-to-leading 
order, while the numerous counterterms of the unnatural parity Lagrangian were neglected. 
The proliferation of such counterterms makes a unique fit to data impossible [4, 7] and one must 
resort to model-dependent assumptions. One possible way of estimating the size of the different 
parameters is to calculate contributions from vector-meson exchange [8], e.g., by employing the 
hidden symmetry formulation of Bando et al. [29, 30]. To this end, an effective Lagrangian of 
unnatural parity and including the nonet of the lowest-lying vector-mesons, V, is constructed. 
Then the vector-mesons are integrated out of the effective theory under the assumption that 
their masses are much larger than the momenta which generates counterterms of order and 
higher. The unknown counterterms of the unnatural parity p^ Lagrangian without vector- 
mesons can thus be written in terms of a few parameters of the vector-meson Lagrangian. The 
coupling constants of the latter can be extracted up to a sign and within certain error bars from 
radiative decay widths of the vector-mesons, such as a; — > 7r°7 and p"*" — > tt'^^. (More recently, 
the parameters of the vector-meson Lagrangian have been constrained by the short-distance 
behavior of QCD [31].) This approach has been applied to the decay process P 77* in [8] 
which amounts to the two-step chain P — > W — > 77* with the virtual vector- mesons being 
treated as infinitely heavy states. Besides giving estimates for the LECs of the higher order 
Lagrangian this procedure reproduces also the experimental P — > 77* slopes [8] and must be 
added to the one- loop contributions to the slopes which are extracted from Eq. (24). These 
are much smaller in magnitude and far away from the experimental data [32, 33, 34]. We will 
discuss this point in more detail in Sec. 5 when we present the numerical results. 

On the other hand, some of the vector-mesons such as the p(770) can be described as bound 
states of two Goldstone bosons [17]. Effective potentials for meson- meson scattering are derived 
from the chiral Lagrangian for the pscudoscalar mesons and iterated within a coupled channel 
BSE which satisfies unitarity constraints for the partial-wave amplitudes. With a small set of 
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Figure 3: Sample rescattering process in the decay P 

parameters a large variety of meson-meson scattering data could be explained up to centcr-of- 
mass energies of 1.2 GeV. In [18] it was shown that the inclusion of the rj' in this framework does 
not spoil the results from [17] for energies below 1.2 GeV as one would expect naively. We adopt 
the approach of [18] here and will not include the vector-mesons explicitly in the theory. They 
will be rather generated from composed states of two pseudoscalar mesons, whereas composed 
states of three pseudoscalar mesons such as the ci;(782) arc beyond the scope of the present 
investigation and will be neglected. The additional parameters in the coupled channel analysis 
are fixed from a fit to the p-wave phase shifts and we do not have to deal with additional coupling 
constants for the vector-mesons which always introduce a theoretical uncertainty. Furthermore, 
our approach is suited to go to higher photon virtualities which presents a problem in the 
resonance saturation scheme, since large momenta in the vector-meson propagators cannot be 
simply dropped when integrating them out of the effective theory. 

Finally, our investigation includes and distinguishes the processes P 1^7'-*^ ^(*)^(*) 
and P VV where the photons can be either on- or off-shell. In the approach 

with explicit vector-mesons, which has been applied so far only to the decays with at least 
one on-shcll photon, the first decay chain is not directly present. It is rather derived from the 
latter one by integrating out the vector-meson coupled to the physical photon. Therefore, for 
the decays with one on-shell photon no clear distinction can be made in this scenario between 
a decay chain with one or two virtual vector-mesons as suggested by complete Vector-Meson 
Dominance. A measurement of the decays with two off-shell photons would help to clarify the 
situation. 

4.1 Bethe-Salpeter equation 

The underlying idea of our approach is as follows. The incoming pseudoscalar meson P can 

decay via a vertex of cither the WZW Lagrangian, Ci^zw, ot the unnatural parity Lagrangian 
at fourth chiral order, C^ct in Eq. (10), directly into one of the following three channels: two 
photons, a photon and two pseudoscalar mesons, or four pseudoscalar mesons. The mesons can 
combine to pairs and rescatter an arbitrary number of times before they eventually couple to 
a photon, see Fig. 3 for illustration. 
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The rescattering process can be described by application of the BSE which generates the 
propagator for two interacting particles. A similar approach has already been successfully 
employed for the hadronic decay modes of r) and rj' [19], where only s-wave amplitudes were 
considered. Here, we extend it to photonic decays and p-waves. In this section, we describe 
the BSE approach for meson-meson scattering in general focusing on the p-wave contributions, 
while in the next section this method will be embedded into the two-photon decays of 7r°, and 
V'- 

In order to describe meson-meson scattering accurately within the coupled channel ap- 
proaches, it is necessary to construct the interaction kernel for the two mesons from the effective 
Lagrangian up to fourth chiral order. In addition to £('^+2) and the Operators O4, O5, On, Ois 
from in Eq. (20), we also include at fourth chiral order the terms 

O23 = (C^) {C,iN) , 02,^ (C^) (C,) (iN) . 

Usually the /3o term is not presented in conventional ChPT, since there is a Cayley-Hamilton 
matrix identity that enables one to remove this term leading to modified coefficients i — 
1, 2, 3, 13, 14, 15, 16 [14], but for our purposes it turns out to be more convenient to include it. 
Hence, we do not make use of the Cayley-Hamilton identity and keep all couplings, in order to 
work with the most general expressions in terms of these parameters. One can then drop one 
of the (3i involved in the Cayley-Hamilton identity at any stage of the calculation. 

Prom the contact interactions of the effective Lagrangian up to fourth chiral order we derive 
the center-of-mass scattering amplitude A{9), where 9 is the c. m. scattering angle. The partial- 
wave expansion for A{9) reads 

2 

= J](2Z + l)AiP^(cose), (41) 
1=0 

where Ai is the effective potential for angular momentum I obtained from the contact interac- 
tions and Pi is the l^^ Legendre polynomial. It is most convenient to work in the isospin basis 
and to characterize the meson-meson states by their total isospin. For the p- waves, e.g., the 
relevant two-particle states have either isospin 0, {KK,riri'), or isospin 1, {117^,1x7], KK,7^rj'). 

For each partial-wave I unitarity imposes a restriction on the (inverse) T-matrix above the 
pertinent thresholds 

ImT.-^ = (42) 

with y/s and being the energy and the three-momentum of the particles in the center-of- 
mass frame of the channel under consideration, respectively. Hence, the imaginary part of T^^ 
is equal to the imaginary piece of the fundamental scalar loop integral Gmm above threshold, 
Eq. (27). 

Following the work of [35, 36] in the baryonic sector, we will adjust the real piece of Gmm by 
introducing a scale dependent constant for each channel in analogy to a subtraction constant 
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Figure 4: Diagrammatic illustration of the Bethe-Salpeter equation for meson-meson rescatter- 
ing. The center-of-mass momentum is denoted by y/s = Qi + (ji = qj + (jf = k + k. 



of a dispersion relation for T^"^. This will compensate the regularization scale dependence of 
Gmfh s-nd bring our results into better agreement with experimental data. In a more general 
way, one could model the real parts by taking any analytic function in s and the baryon and 
meson masses. This option has been successfully applied for the case of SU (2) baryon ChPT 
in [37], but is not necessary in the present work, as we can reproduce the experimental phase 
shifts for s- and p-wave meson-meson scattering already very accurately with the first method. 
The inverse of the amplitude can be decomposed into real and imaginary parts ^ 

T-^ = T-^ + G (43) 

with 

Gmm = GmfhilJ') + amm(A^), (44) 

where and Re[G'] give the real part and Im[G'] gives the imaginary part of T^^ as required 
by unitarity, Eq. (42). The scale dependences of G and a on /i cancel each other and we will 
choose the constant a to depend on the angular momentum /. Inverting (43) yields 

T^[l+r-G]-^r (45) 

which is understood to be a matrix equation that couples the different channels. The matrix 
G is diagonal and includes the expressions for the loop integrals in each channel. Expanding 
expression (45) 

T^r-T-G-T... (46) 
and matching the first term in the expansion to our tree level amplitude for each partial-wave 

Ti = Ai (47) 

our final expression for the T matrix reads 

T=[l + A- G]-^ A (48) 

which amounts to a summation of a bubble chain in the s-channel. This is equivalent to a 
Bethe-Salpeter equation with A as potential, see Fig. 4. 

The unknown LECs of the chiral Lagrangian must be fitted to experimental data. This has 
partially been accomplished in [18], where after applying the same approach as in the present 
investigation we constrained the LECs of the Lagrangian up to fourth chiral order by comparing 
the results with the experimental s-wave phase shifts of meson-meson scattering. Agreement 

^For brevity we suppress the subscript I. 
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Figure 5: Fit to the experimental phase shifts ([38, 39]) for tttt — > tttt in the I — J — 1 channel. 

was achieved in the isospin J = 0, | channels up to 1.2 GeV and in the isospin / = |, 2 channels 
up to 1.5 GeV. However, variations of some of the parameters which have been set to zero for 
simplicity do not yield any significant effect for the s-wave phase shifts. In fact, some of them 
could be constrained from the hadronic decays of the rj and rj' [19]. A good fit to the decays 
[19] and the s-wave scattering data [18] is given by 

^« = = 0, 

(3^ = 0.56 X 10-3, = -0.3 X 10-3, (49) 

/5f = 1.4 X 10-3, = 0.06 X 10-3 

with v^^ = \P — V^'^a^'* — 3^2^^ and all the remaining parameters being zero. It is not trivial 
that with such a small number of parameters we have been able to explain a variety of data 
within the approach. 

In the present investigation, we are particularly interested in the p-wave amplitudes which 
will be the only contributions to the photonic decays. Our fit to the scattering data from [38] 
and [39] is shown in Fig. 5. We are able to use the same set of parameters as for the s-wave 
scattering in [18] with only one non-vanishing scale dependent constant 

am^m-. = -4.37 X 10-2 (50) 

at /X = 1.0 GeV and we take ji = 1.0 GeV in all channels. It is also important to note, that 
the parameter choice in Eq. (49) is not unique, since variations in one of the parameters may 
be compensated by the other ones. Nevertheless, we prefer to work with this choice; it is 
capable of describing the experimental phase shifts both for s- and p-waves with a small set 
of parameters and motivated by the assumption that most of the OZI- violating and l/N^- 
suppressed parameters are not important, although /3g°'* and v'^\v2^ in V2\v2'^ have small but 
non- vanishing values. 

4.2 Meson-meson rescattering in tt^, r/, r]' — > 77 

The meson-meson rescattering processes from the previous section can be employed as effective 
interaction kernels for two mesons in the decay process, see Fig. 6. In order to perform the re- 
maining loop integrations in Fig. 6 which are not included in the effective two-meson interaction 
kernel iT, we rewrite the partial-wave decomposition of the T matrix of the meson-meson scat- 
tering process as follows [18] . The general expression for T depends only on scalar combinations 
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Figure 6: Set of meson-meson rescattering processes in the decay P ^ included in our 

approach. 



of the momenta which can be expressed in terms of the Mandelstam variables. The Mandel- 
stam invariants s, t and u are defined as the center-of-mass energy squared s = (g^ + qif = p'^, 
the momentum transfer squared t = {qi — qfY and the crossed momentum transfer squared 
u — {qi — qfY, where we used the notation of the previous section, see Fig. 4. The constraint 
s + t + = + + + ^1 = + + + allows one to neglect the combination 
t-\-u m favor oit — u and the scalar amplitude can be written as T(s, t — u). Since we restrict 
ourselves to the effective Lagrangian up to fourth chiral order, the partial-wave decomposition 
of T — which is the expansion of T in t — u — is given by 



T — Ylii "^iJi — TgJs + TpJp -\- TdJd, 



(51) 



where the partial-wave operator Ji is a polynomial of degree I in t — u and proportional to the 
Legendre polynomial Pi in the cosine of the scattering angle. The J; can be written as 



J. 



Jp = hi^.^q^q'} = — — + 



t-u {q- - qi){qf - Qf) 



hixvqtqi hpaqfqf 



J,. =Jt - 



4s 



(52) 



with 



(53) 



It is now straightforward to perform the remaining loop calculations in Fig. 6, since the TJ 
depend only on the Mandelstam variable s, i.e. the momentum squared of one of the photons. 
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Hence, they are independent of the loop momenta and can be factored out of the loop integral. 
The remaining pieces Ji can be expressed in terms of the loop momenta and the loop integrations 
can be performed. In this way, the T matrix is treated as an effective vertex, which summarizes 
the infinite chain of meson-meson rescattering processes. One checks explicitly that s- and 
(i-waves drop out, whereas only the t — u piece in Jp contributes. Only the channels (tttt, KK) 
contribute so that = g^, and hence the limit s for on-shell photons does not yield a 
divergence in Jp. The results for the diagrams of Fig. 6 read 



X 



+i{ml; k'^)[fir^''^\k'^)i{ml; k'^) + f^''^''^\k'^)i{ml; k'^)] 
for the one vector-meson exchange and 

X i{ml- k^)[f^''-''^\k^)i{ml- k') + f^-''''\k')i{ml; k')] 



(54) 



(55) 



X I{ml;k'^)[fi'-''^\k'^)Iiml-k'^)+t^^^^ 

for the two vector-meson exchange, where Tp'^^'^^ denotes the p-wave amplitude in Eq. (51) 
for channel a scattering into channel b. The symbol E' indicates summation over the meson 
pairs TT+TT^, K^K~ and K^K^. Note that the exchange of two vector mesons arises due to the 
five-meson vertex in the piece (E^) of the WZW action, Eq. (6). The coefficients in Eq. (54) 
are given by 

1 

2' 



(1 CC),7r± _ (1 CC),K± _ (1 CC),K^K° 
- %o - %o 
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+ 4/3, 



5,18 



The coefficients ry(^ '^'-') (.^sg of two coupled channels are symmetric under a b. The 
non- vanishing ones read 



(2CC),iv'^,K^ _ (2CC),-K^,K°K° _ 3 
77r0 — 77r0 — ^ ) 



(2CC),-K^,K^ _ (2CC),Tr^,K°K° _ 1 (2 CC),K'^ ,K°K" _ 



Iri' 



(2CC),%'^,K^ 



(2CC),Tr^,K°K° 



'2% 



(2CC),K^,K°K° 



4 ' 

2 , 



m 



K 



ml) 



.,(2) 



+ 4/3. 



'5,18 



(57) 
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The loop integral / in Eqs. (54,55) is given by 

where the integral / is defined in Eq. (25). Here we made use of the freedom to take arbitrary 
values for the analytic pieces of the integrals which corresponds to a specific choice of coun- 
terterm contributions. To be more precise, we have kept in the pion loops a term of the type 

(Cjr — while neglecting all other analytic portions. Alternatively, one could have 

altered the regularization scale // of the integrals relevant for the term, but we preferred to 
take explicit analytic pieces, while keeping ^—1 GeV fixed in all channels, which is similar to 
adding a subtraction constant as in Eq. (44). As we will see in the next section, the inclusion 
of such analytic portions which are beyond the accuracy of the one-loop calculation yields an 
improved fit to the experimental data and accounts for dynamical effects of higher chiral order. 

Clearly, we are missing further unitarity corrections beyond the one-loop calculation. How- 
ever, from the discussion below it will become clear that the set of diagrams included in our 
model is consistent with available data. 



5 Numerical results 

In this section, we will discuss the numerical results of our calculation. In order to compare 
our results with experimental data, we include the counterterm contributions from Eq. (34). 
At the one-loop level, it was not necessary to take them into account, as the one- loop formulae 
could already be brought to agreement with the experimental decay widths by fitting either 
Fy^i and v^2^ or v^'' and one of the contact terms w^, w^'' ■ The inclusion of the coupled 
channel analysis, on the other hand, leads to changes in the decay amplitude for two on-shell 
photons which must be compensated by counterterms that are proportional to the quark mass 
matrix, i.e. those with the coefficients wf\ wf\ w^^ and w^^\ It turns out that the inclusion 
of counterterms of sixth chiral order as discussed in Section 3.1 is sufficient to compensate these 
additional contributions. Furthermore, A;^ dependent counterterms {wf'\ w^"*) arc needed to 
obtain agreement with data for space-like photons with squared four-momenta A;^ < 0, since 
the non-analytic contributions from the one-loop diagrams are too small to account for the 
behavior of the transition form factor in the space-like region, where furthermore contributions 
from the BSE are almost negligible. The two possible terms from the sixth order Lagrangian 
C^ct 1 Eq- (33), are sufficient to bring our results into better agreement with the data for photon 
virtualities up to k"^ pa —0.8 GeV^. The values of the counterterms are (in units of GeV~^) 

wf'^ = -5.5 X 10-^ w^2^ = -1.5 X 10-^ wf = 1.31 x 10"^ , . 

ivP = -1.4 X 10-^ lyf ^ = -0.47 x 10-^ w^^^ = -0.16 x 10"^ ^ ^ 

The parameters w^^ and Wg^^ enter only in the combination |(m|; -|- 2ml)w^^^ + A(2m\- + 
m^)wg^'', see Eq. (35), hence their values cannot be fixed separately. Since the Wg^^ term is 
1/A^c suppressed, we choose w^^^ — wf^ /?>. 

The results of our model are compared in Fig. 7 with the transition form factors for the 
decay of the 7r°, and r]' into one on-shell and one off-shell photon. The transition form factor 
Tp{k'^\k''^) is defined as 

^(P ^ ^*^(*)) ^ k.e^k'Js^^'^'^^Tpie; k") (60) 
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Figure 7: Results of the one-loop calculation (dotted), including the diagrams Figs. 6a+b 
(dashed), full calculation (sohd). Data: sohd squares [40], open squares [41], diamond [26], 
open circles [42] , sohd circles [43] . 



and we plot the quantity ^\J^p(k'^;k''^ = 0)P which yields the width of the decay into two 
physical photons at k'^ = 0. 

The overall good agreement with the data indicates that our model is capturing the impor- 
tant physics. In our fit^ we have set the non-anomalous terms of unnatural parity at fourth 
chiral order to zero, — W2'^ — w^^ — 0, since a good fit to the data can already be 
achieved without them, while keeping F^i fixed at F^i — 1.1 F^. Furthermore, we set the 
mixing parameter v^'^ to zero which is consistent with a previous coupled channel calculation 
on hadronic decays of the rj and rj' [19]. It should be emphasized that for a non-perturbative 
coupled channel approach the values of the coupling constants do not necessarily coincide with 
those from a perturbative loop-expansion, which was already pointed out in [19]. It is therefore 
not surprising that v^^ and the w^^^ differ in both schemes. The inclusion of the diagrams in 
Figs. 6a, 6b which contain only one coupled channel and thus are associated with the exchange 
of only one vector-meson already yields the crucial structure of the curves (dashed line). The 
exchange of two vector-mesons (Fig. 6c) is merely a small correction which could even be com- 
pensated by a change of the parameters of the 0{p^) counterterms. For completeness we have 
also shown the numerical results from the one-loop calculation which have been supplemented 
by the A;^ dependent counterterms at sixth chiral order, in order to be in better agreement for 
k'^ < 0. Due to this procedure the k"^ dependent terms differ from those in the coupled channel 
approach. Nevertheless, with this choice of parameters it is not possible to reproduce the sharp 
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Figure 8: Dependence of JF,, and J^rf on the mixing parameter v^^: = (solid), — 
0.5F2/4 (dashed), = 1-15^2/4 (dotted). 
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Figure 9: Dependence of JT^/ on the LEG w^^^: w^^^ = -3.0 x lO"^ (dashed), w^^^ = (sohd), 
w^^^ = 3.0 X 10-=^ (dotted). 



increase of the form factor for time-hke photons, > 0, at the one-loop level, and we miss in 
any case the resonance structure in the f]' decay. 

The dependence of our results on the rj-rj' mixing parameter v!^^ is depicted in Fig. 8. We 
have plotted the curves for different vahics of v^2 \ which indicate that a value for ■52^'' around 
zero is in better agreement with the data. The results for 7r° —>■ 7*7 are, of course, not affected 
by changing {ig^^ 

In Fig. 9 we show our results for different values of the LEG When keeping v^^ — 
only the process r)' — > 7*7 is altered. From the plot it can be seen that ^3^^ has significant 
influence on the peak structure at k'^ ~ 0.6 GeV^ and we can exclude values tOg^^ > 2.0 x 10~^. 

Another quantity of interest is the slope parameter of the transition form factor. It is defined 
as the logarithmic derivative of the form factor jFp at /c^ = 

(61) 

In Table 1 we compare our values with simple pole fits to the form factors. The value for 
7r° is the world average given in [26]. For r] and rj' we cite three measurements, the Lepton-G 



bp= :77^1n^p(A;';A;'2 = 0) 
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Experiment 


NLO 


1 CC 


full 




1.76 ±0.22 [26] 


1.19 


1.98 


1.95 


7] 


1 9 ± 4 [421 

J. • t_/ 1 \J m 1 -—J 1 

2.0 ±0.5 [41] 
1.42 ±0.21 [40] 


0.73 


1.57 


1.58 


T)' 


1.7 ± 0.8 [43] 
1.38 ±0.23 [41] 
1.59 ±0.18 [40] 


0.85 


1.78 


1.79 



Table 1: Slope parameters bp in units of GeV~^ derived from pole fits to the experimental data, 
the next-to- leading order calculation, including one coupled channel (diagrams Fig. 6a±b), and 
the full calculation. 

Collaboration [42, 43], the TPC/27 Collaboration [41] and the CELLO Collaboration [40]. From 
a pole fit a characteristic mass Ap is extracted and related to the slope parameter via bp = 1/ Ap. 
The obtained values arc given by A^ = (0.72 ± 0.09) GeV [42], A^ = (0.70 ± 0.08) GeV [41], 
A^ = (0.84 ± 0.06) GeV [40] and A^/ = (0.77 ± 0.18) GeV [43], A^, = (0.85 ± 0.07) GeV [41], 
Aj]' — (0.79 ± 0.04) GeV [40]. With one exception we achieve good agreement within error bars 
with the pole fits to the form factors. We point out that in the determination of the given 
errors for bp it has not been taken into account that in principle a different ansatz for the form 
factor can lead to different values of the slope parameter. In our model a significant part of the 
slope parameters is induced by the meson-meson rescattering processes. Again the exchange of 
two vector-mesons plays a negligible role. 

Within our model, contributions from composed states of three pscudoscalar mesons which 
would correspond to the c<j(782) have been neglected. However, the effects of the u;(782) should 
be similar to those of the p(770) due to their almost equal masses. It may well be, that by 
keeping some of the analytic portions of the loop integrals in Eqs. (54,55) which amounts to 
a particular choice of higher order counterterms, we took into account the dynamical effects 
of the a; (782). Within our model, the contributions from the vector-mesons should thus be 
regarded to be a combination of both the a;(782) and the p(770). 

The contributions of the two vector- meson exchange diagrams, on the other hand, are 
almost negligible, as can be seen from Fig. 7. This is in sharp contradistinction to the complete 
Vector- Meson Dominance picture where the photons can only couple via the exchange of vector- 
mesons. Of course, when coupling to an on-shell photon the pertinent vector-meson propagator 
reduces to a vertex, but with both photons off-shell, the two approaches should yield different 
predictions. In Fig. 10 wc show our predictions for two off-shell photons for 11^,1],')]' decays 
with A;^ = k''^ as well as the rj' transition form factor in the (/c^, /c'^)-plane. As in the case with 
one off-shell photon the result of a calculation with one coupled channel is almost identical 
to the full calculation. A measurement of the transition form factors for two off-shell photons 
could serve as a check of our model and help to clarify, whether the exchange of only one or 
rather two vector-mesons occurs in these decays. For one and two off-shcll photons with larger 
space-like momenta the pion transition form factor has been calculated within an instanton 
model of QCD [44] . For small negative photon virtuahties the results are quite similar to ours. 
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Figure 10: Prediction of the transition form factors for two off-shell photons with A;^ = A;'^ 
resulting from the full calculation (solid) and from a calculation without the contribution from 
Fig. 6c with two coupled channels (dashed, almost on top of solid line). For Vj 7*7* we also 
show the transition form factor from the full calculation in the (A;^, A;'^)-plane. 



23 



6 Conclusions 



In the present work, we have investigated the two- photon decays of tt*^, 77 and r( within a 
chiral unitary framework. To this end, we have supplemented the one-loop calculation of chiral 
perturbation theory by a Bethc-Salpctcr approach which satisfies unitarity constraints and 
generates vector-mesons from composed states of two pseudoscalar mesons. 

While the one-loop result is sufficient to achieve agreement with the decay widths P ^ 77 
(P = 77^,77,77'), the vector- meson exchange is crucial to reproduce the sharp increase of the 
transition form factor for time-like photons in the decays P 77*. Our method reproduces 
also the resonance structure in the transition form factor of the decay 77' 77* at photon 
virtualities around /c^ ~ 0.6 GeV^. 

Furthermore, our approach distinguishes between single and double vector-meson exchange, 
the latter one being the only contribution in the complete Vector-Meson Dominance picture. 
Our study suggests that for the decays with off-shell photons the exchange of one vector- 
meson is the dominant contribution, whereas the two vector-meson exchange is almost negligible 
which is in contradistinction to complete Vector-Meson Dominance. The available data on 
electromagnetic transition form factors is restricted to decays with exactly one off-shell photon. 
When coupling to an on-shell photon, the pertinent vector-meson propagator reduces to a 
vertex so that no clear distinction can be made between the two approaches. However, an 
experiment with two off-shell photons should be able to clarify, whether one or two vector- 
meson exchange occurs. The question whether double Vector-Meson Dominance holds is also 
an important issue for kaon decays and the anomalous magnetic moment of the muon. We have 
presented predictions for the decays with two off-shell photons, and again the exchange of two 
vector-mesons plays a minor role. 

Within our model, contributions from composed states of three pseudoscalar mesons which 
would correspond to the cj(782) have been neglected. However, the effects of the cj(782) should 
be similar to those of the p(770). It may well be, that by keeping some of the analytic portions 
of the loop integrals in our coupled channel analysis which amounts to a particular choice of 
higher order countcrtcrms, we took into account the dynamical effects of the cl^(782). Within our 
model, the contributions from the vector-mesons should thus be regarded to be a combination 
of both the tu(782) and the p(770). 

The present approach can be applied in a straightforward manner to the anomalous decays 
77, 77' — > 7r'^7r~7 and the radiative decays 77, 77' 7r'^77 which will be investigated in future 
studies. 
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A Counterterms of order 0{p^) 

In this section, we discuss the relevant counterterms of order 0{p^). For the covariant derivative 
and the field strength tensors we use the following notation 

D^U = d^U-tr,U + tm^, ^^^^ 
Rfjiv dijTi, di/fn ^[^/i) ^i/]) -^/ii/ dijlu di/l^ i^fjii 

where = f/i + Oju, — — a^^ and dfj, — (a^) . The replacement of the original 

field strength tensors R^^, L^i, by the new quantities R^^, L^,^ leads to additional counterterms 
which involve the derivative of the singlet component of the axial- vector field, {d^a„), however, 
such counterterms do not contribute to the processes discussed in the present work and are 
neglected. 

In the SU (3) framework there are numerous terms containing six derivatives which according 
to [24] can be reduced to only one contribution by utilizing methods such as partial integration, 
equation of motion, epsilon relations and Bianchi identities. In our notation it is given by 

i Wie^"'^" {{2U[D^Rx^U + 2D^Lx^ + WRx^D^U + D^WRx^U + UW^ULx^W + Lx^^D^WU) 

(A.2) 

where we made use of the abbreviations 

H^u = U^R^.uU + L^,, C^ = WD^U. , . 

M^Wx + X^U, N^Wx-X^U. ^ ■ ' 

For the decays into two photons Eq. (A.2) yields the contribution (in the differential form 
notation of [5]) 

-yj(^) }^ (^(f)dvDdv). (A.4) 

From the extension to the U (3) framework there arise several more terms which reduce after 
applying the above mentioned methods to the relevant contribution 

wi^^ J Vo (dvndv). (A.5) 

Furthermore there are terms of order 0{p^) that contain the chiral symmetry breaking 
object X = 2BM. with Jv[ being the quark mass matrix 

The first two terms arise from conventional SU (3) theory whereas the last two are due to the 
extension to U{?>). In the differential form notation the relevant contributions to the decay 
P — > are given by the contact terms 

4S = '^{<f^xdvdv) + '^{<l^x){dvdv) 

+ w'i'' ^ 7^0 (x dv dv) + w^Q^ ^ 7^0 (x) {dv dv). 
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Combining all counterterms of order 0{p^) we obtain 

£g) = -wf ) 1^ (0 dv □ dv) + w^'^ y r^o {dv □ rf^;) + ■ (A.8) 
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